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As the comment of Tiwari [1] reflects typical and prob- 
ably common misunderstandings, we feel it valuable to 
elaborate further on these delicate and important issues: 

(i) In the gauge coupling vpj^A^tp, the gauge field A^ 
play a dual role: it provides a physical coupling to the 
Dirac field vb, as well as a gauge freedom to compensate 
for the phase freedom of i\>. Our idea of solving the gauge- 
invariance problem is to decompose this dual role by seek- 
ing, in any gauge, a unique separation A = A pure + v4 p h ys , 
with Apurc a pure-gauge term transforming in the same 
manner as does the full A and always giving null B, and 
^.phys a physical term transforming in the same manner 
as does the electric field E. Namely, A p h ys is gauge in- 
variant/covariant in electrodynamics/Yang-Mills theory, 
while A puro has the same gauge freedom as A and can 
be used instead of A to construct a covariant derivative 
-Dpurc = V - iA puic . The separation A = A pmc + A phys , 
and the definitions for A p h ys and A puro , are by no means 
specifying a gauge or restricting the gauge freedom of 
A. E.g., in electrodynamics, V • ^4 p h ys = is the defini- 
tion of j4 p hys in any gauge, it is not to be confused with 
the gauge condition V • A = 0. Nevertheless, the sepa- 
ration A = Apurc + Aphys does become the simplest in a 
unique physical gauge in which A pure = and ^4 p h ys = A. 
Coulomb gauge V- A = is the physical gauge for electro- 
dynamics. But for Yang-Mills theory the physical gauge 
is [A, E] = which makes the gauge-dependent gluon 
color charge vanish [2| , while V ■ A = is no longer priv- 
ileged. This may elucidate why quantization of Yang- 
Mills theory in the Coulomb gauge is not very illuminat- 
ing. 

(ii) Our solution achieves the best possible reconcil- 
iation of Lorentz covariance and gauge invariance for 
the momentum and angular momentum (AM) of cou- 
pled quarks and gluons. To appreciate this point, we 
first note that, in the instant form, six Poincare gener- 
ators can be interaction- free (or good): P = P q + P g , 
J = J q + J g , and four generators involve interactions: 
H = H q + Hg + H int , K = Kg+Kg+Kint. Here q, g, and 
int denote the quark, gluon, and quark-gluon interacting 
parts, respectively. In gauge theories, it should be further 
noted that not only the interaction-involving Lorentz 
transformations, but also the interaction-free ones, are 



troublesome. This is because an interaction-free opera- 
tor, like J d 3 xyj^xx jVV>, is often gauge-dependent, while 
a gauge-invariant construction, like J d 3 xijj^xx iDi/j, of- 
ten involves interaction. In fact, the essence of our con- 
tribution is that the good generators P and J are in- 
deed constructed to be both interaction-free and gauge 
independent, so that each part in P and J transforms 
properly under spatial translations and rotations. But 
since K involve interaction intrinsically, only the total 
J transform properly under boost: [K l , J J ] = ieijkK k , 
while [K l , J| fl ] = ieijkK* can not hold (otherwise we 
would have [K l , J° q + Jj] = iey* [K* + Kg] , which contra- 
dicts [K l , J J ] = ieijkK h ). To determine the boost prop- 
erties of J q ,g, we need to carry out a canonical quantiza- 
tion (preferably in the physical gauge so that the pure- 
gauge term vanishes), then compute the commutators 
[K l ,Ji ]. This task is difficult but unavoidable for ev- 
eryone. In comparison, the light-cone formalism renders 
boost along the third axis interaction-free, but makes the 
description of AM rather troublesome, because two AM 
components (namely, the rotation generators along the 
x, y axes) involve interactions. All in all, the separation 
of P or J into quark and gluon parts can never be fully 
covariant under all Lorentz transformations, no matter 
how one deals with the gauge field. In fact, this property 
is an intrinsic complication for any interacting system, 
not merely gauge interactions. By decomposing A into 
physical and pure-gauge terms, we add absolutely no ex- 
tra complication in this regard. 

(hi) It is an illusion that the spin and orbital AM of 
light as in the optic measurements or multipole-radiation 
analysis can be straightforwardly interpreted in classical 
electrodynamics, without applying the approach which 
we follow. As long as one discusses spin and orbital AM 
separately, the gauge-invariance problem is sharply en- 
countered, no matter one adopts a classical or quantum 
language. This problem is just absent if one merely con- 
cerns about the integrated total AM. Without defining 
the spin and orbital AM gauge-invariantly, it is impos- 
sible to interpret the optic measurements which intend 
to manipulate spin and orbital AM separately, or the 
multipole-radiation analysis which employs the notion of 
spin-orbital coupling. We should also clarify that the so- 
lution we provide applies to both classical and quantum 



dynamics. It is our gauge-invariant construction for the 
spin and orbital AM that justifies the spin and orbital 
values or quantum numbers assigned to the electromag- 
netic field or photon in any experiment. And in turn, the 
consistency of our theory with experimental results sup- 
ports the validity and correctness of our gauge-invariant 
expressions. 

(iv) To see the physical content of E x B, note that 
J = J d 3 xx x (E x B) gives the total AM of a free gauge 
field, including both spin and orbital parts. But if P(x) 
is the momentum density, then J d 3 xxx P(x) is the stan- 
dard form of orbital AM. This implies that Ex B is not a 
pure mechanical momentum, it must include a spin flow 
so as to give the total J by an apparent orbital form. By 
writing J = f d 3 xE x A phys + J d 3 xx x (E^A^A, the 

correct momentum density is read out to be E l VA l h 
From ExB = E l VA phys - V 4 (£M phys ) + (V • E)A phys , 

we see that the difference between E x B and E l VA l h 
is a surface term only if the (in general charged, here) 
Dirac field is absent so that V • E — 0. For an inter- 
acting system, the difference is substantial. Similar sit- 
uation occurs for the AM: x X (E X B) = E X A p h ys + 
x x (E*VA phys ) + V l [E l (A phys x x)] - (V • E)A phys x x, 
hence f d 3 xx x (E x B) may differ drastically from 
/ d 3 xE x Aphy S + / d 3 xx x (E l VA phys ) in the presence of 
interaction. For the nucleon, then, E x B and E l VA 1 ^ 



may give totally different views of how much nucleon 
momentum is carried by gluons. E x B is a compo- 
nent of the symmetric energy-momentum tensor, while 
E*VA l h belongs to the (non-symmetric) canonical one. 
These two tensors differ by a total-derivative spin term. 
This term does not matter for the integrated total energy 
or momentum, so there is no curiosity in that a photon 
with arbitrary polarization and orbital AM can have the 
same energy. However, the spin term does alter the den- 
sity and symmetry of the energy-momentum tensor, for 
which we do need a concrete density expression for the 
purpose of coupling to gravity. Thus, were one to adopt 
Einstein's gravitational equation which requires a sym- 
metric tensor, then the "momentum" that contributes to 
gravity is not the mechanical momentum [3| . 

In closing, we remark that Ref. [2| does solve the long- 
standing gauge-invariance problem of spin and orbital 
AM, in the sense that the expressions satisfy all theoreti- 
cal requirements and agree with all experimental results. 
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